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Noise correlations are studied for systems of hard-core bosons in one-dimensional lattices. We use 
an exact numerical approach based on the Bose-Fermi mapping and properties of Slater determi- 
nants. We focus on the scaling of the noise correlations with system size in superfluid and insulating 
phases, which are generated in the homogeneous lattice, with period-two superlattices and with 
uniformly distributed random diagonal disorder. For the superfluid phases, the leading contribution 
is shown to exhibit a density- independent scaling proportional to the system size, while the first 
subleading term exhibits a density-dependent power-law exponent. 
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I. INTRODUCTION 

In the past decade, ultracold quantum gases have 
gained considerable attention due to the unique control 
achieved experimentally for manipulating such systems. 
This has enabled experimentalists to explore the richness 
and complexity of strong correlations and reduced dimen- 
sionalities and to even simulate model Hamiltonians used 
to understand complicated materials [l| . For the goal of 
studying quantum systems in quasi-one-dimensional ge- 
ometries, remarkable experimental examples include the 
realization of quantum gases in very anisotropic traps 
[U, HI and loading Bose-Einstein condensates in deep two- 
dimensional optical lattices [l-Q and in atom chips [Tol — 
EH- 

One dimension hosts a variety of models that can be 
exactly solved analytically and as such are of very much 
interest to both theorists and experimentalists. Remark- 
ably, the high degree of tunability and isolation achieved 
in ultracold gases experiments has permitted the realiza- 
tion of various such models. An example of particular 
relevance to the work presented here was the realization 
of a gas of impenetrable bosons (hard-core bosons), also 
called a Tonks-Girardeau gas, in the presence Q and ab- 
sence d, @ of a lattice along the one-dimensional gas. 
The problem of indistinguishable impenetrable bosons in 
one dimension was first analyzed by Girardeau [l3j , who 
noticed that its thermodynamic properties could be eas- 
ily computed by mapping such a problem to that of in- 
distinguishable noninteracting spinlcss fcrmions. In the 
presence of a lattice, the hard-core boson problem (see 
below) can be mapped to a special case of the XY spin- 
1/2 chain introduced by Lieb et al. whose thermody- 
namic and local properties can also be solved by mapping 
it to a noninteracting spinless fermion lattice model. 

The calculation of the off-diagonal correlations, such as 
the one-particle correlations, is a more challenging task. 
In the homogeneous case, this has been done in various 
works and using variou s ap proaches for both continuous 
and lattice systems [15l - [20( ] . It should be noted, however, 
that the experimental realization of these model Hamilto- 
nians requires a trap for containing the gas. This means 



that such experimental systems are in general inhomoge- 
neous and their description requires one to take into ac- 
count the presence of the trapping potential, which is to 
a good approximation harmonic. Studies of one-particle 
correlations of harmonically trapped Tonks-Girardeau 
gases have been performed in a series of more recent 
works [M^27|. 

One-particle correlations can be probed in experiments 
by means of time-of-flight measurements, in which the 
confining potentials arc turned off and, in the absence 
of interactions during the expansion, the initial momen- 
tum distribution of the trapped gas is mapped onto the 
density distribution of the system after a long expansion 
time. The latter density distribution is then measured by 
taking a picture of the gas after expansion. How the scal- 
ing of the one-particle correlations in the trapped system 
is reflected in the momentum distribution, which is the 
diagonal part of the Fourier transform of the one-particle 
density matrix, was also discussed in several works men- 
tioned above ■ 

Remarkably, it was also proposed that higher order cor- 
relations can be measured after time of flig ht by analyz- 
ing the atomic shot noise in the images [28[ ■ These noise 
correlations are experimentally associated with Hanbury- 
Brown-Twiss interferometry, which allow one to measure 
the density-density correlations in the spatial images. Af- 
ter long expansion times, under the usual assumption of 
absence of interactions during the expansion, noise corre- 
lations reflect the momentum space density-density cor- 
relations in the t rapp ed system. Shortly after the the- 
oretical proposal [231, noise correlations were measured 
in experiments with bosons in three-dimensional optical 
lattices [2!| and with attractive fermions [3fjj . 

Our goal in this paper is to explore the scaling of the 
noise correlations in various ground-state phases of one- 
dimensional hard-core-boson-lattice systems. We con- 
sider the homogeneous case, systems with an additional 
period-two superlattice potential, and disordered sys- 
tems with a uniform random distribution of local po- 
tentials. We implement an exact numerical approach to 
compute the noise correlations, which follows after the 
hard-core-boson-latticc model is mapped onto a noninter- 
acting spinlcss fermion model by means of the Holstein- 
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Primakoff transformation [3l| and the Jordan- Wigner 
transformation [32"j . This approach is an extension of the 
method developed by one of the authors (in collaboration 
with A. Muramatsu) [111 [26[ for the exact calculation 
of the one-particle density matrix of hard-core-boson- 
lattice systems using properties of Slater determinants. 
We should note that earlier studies of noise correlations 
in hard-core-boson-lattice models followed an alternative 
numerical formulation based on Wick's theorem [33j - l35j . 
However, the lattice sizes accessible within that approach 
were too small to enable a systematic study of the scaling 
of the noise correlations with system size. 

There are three ground-state phases on which we will 
focus our present study, which are the superfluid phase, 
the Mott-insulating or charge-density-wave phase, and 
the Anderson-glass phase. Those phases can be obtained 
in the various background potentials mentioned before. 
In the superfluid phase, we show that the leading con- 
tribution to the noise correlation peaks scales linearly 
with the size of the system, independently of the density 
and of the absence or presence of a supcrlattice poten- 
tial, while the first sublcading term does depend on both. 
As expected, for the Mott and Anderson-glass phases, 
which are both insulating, the scaling of the peaks shows 
an asymptotic value that depends on the density and 
strength of the background potential but that is inde- 
pendent of the system size. The leading-order results 
are consistent with the behavior of the zero-momentum 
peak of the momentum distribution, which scales with 
the square root of the system size in the superfluid phases 
[l5j while it saturates in the insulating (3(| H3] and dis- 
ordered [38[ phases. The latter behavior is a result of the 
short-range correlations present in the insulating phases. 

This presentation is organized as follows. In Sec. [Til 
we describe the models and introduce the exact numeri- 
cal approach. In Sees. IIIIl IIV1 and El we study the noise 
correlations in the homogeneous case, in period-two su- 
pcrlatticcs, and in disordered systems, respectively. A 
comparison between the noise correlations in all those 
systems is also presented in Sec.[V] Finally, Sec. |VI] sum- 
marizes our results. 



II. EXACT APPROACH 
A. Hamiltonian and relevant quantities 



[b i , Sj] = Sij still hold for all sites, additional on-site con- 
straints apply to the creation and annihilation operators 

&r=^ , (2) 

which preclude multiple occupancy of the lattice sites. 
Note that Eq. $2$ is only valid when applied to a string 
of bosonic operators in normal order [33j , as will be ex- 
plained below. 

The hard-core-boson Hamiltonian can be mapped onto 
the exactly solvable noninteracting fcrmion Hamiltonian 
by means of Bose-Fermi mapping, which follows in two 
steps. The first step is given by the correspondence be- 
tween hard-core bosons and spin- 1/2 systems through the 
Holstein-Primakoff transformation [3l| 

<t+ = b\ y/l-bl'bi, gT = \j\-b% 

&i = bfa-l, (3) 

where af 1 are the spin raising and lowering operators and 
of is z-component Pauli matrix for spin-1/2 systems. A 
straighforward analysis reveals that can be directly 

replaced by (a~ ) if and only if the hard-core boson cre- 
ation and annihilation operators are arranged in normal 
order; that is, all creation operators must be placed to 
the left of the annihilation operators before the mapping. 
The root of this difference between hard-core-boson and 
spin-1/2 systems lies in the fact that despite the sup- 
pressed multiply-occupied states, virtual states of mul- 
tiple occupancy can occur in the infinite U limit of the 
Bosc-Hubbard model and they need to be properly taken 
into account for a correct calculation of bosonic correla- 
tions [33j . As mentioned above, in general Eq. (J2) docs 
not apply, for example, for a bosonic system (indepen- 
dently of the value of U): (0|66&t&t|o) = <1[6 6+|l) = 2, 
and a direct replacement of the hard-core-boson opera- 
tors by spin operators would lead to a strictly zero ex- 
pectation value. To avoid this problem, the proper recipe 
is to normal order strings of hard-core-boson operators 
using the bosonic commutation relations before making 
the replacement St(6 i ) — > crf(a^). 

In the second step, the spin-1/2 Hamiltonian can be 
mapped onto a noninteracting fcrmion Hamiltonian by 
means of Jordan- Wigner transformation [32| . 



In the hard-core limit of the Bose-Hubbard model, the + ~t TT -ivpf »- TT inPf " 

one-dimensional Hamiltonian can be written as a i h 11 e p ■> ®i = 11 e ^ 13 hi 

13 = 1 p=l 

Hhcb = -tJ2 (Wh+i + H.c) + W &* = f}f. - I, (4) 

i i 2 



where t represents the hopping parameter and {V{\ a set 
of on-site potentials. The hard-core boson creation and 
annihilation operators at site i are denoted by b\ and b i , 
respectively, and hi = 6^6^ denotes the occupation oper- 
ator of site i. While the bosonic commutation relations 



with fl and f t being the creation and annihilation oper- 
ators for spinless fermions, respectively. 

The noninteracting fermions share the exact same form 
of the Hamiltonian as the hard-core bosons up to a 
boundary term that for periodic systems depends on 
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whether the total number of bosons N in the system is 
even or odd 14011: 



(5) 



where hj = fj f i is the fermionic occupation operator 
of site i. This mapping shows that all thermodynamic 
properties and real space density-density correlations of 
hard-core bosons are identical to those of a system of 
noninter acting fermions. This is of course not true for 
the off-diagonal correlation functions. 

In order to compute the one-particle correlations, one 
can follow the approach described in Refs. |25|,[26j]. (Note 
that in those studies the hard-core boson and spin-1/2 
operators were used indistinctively but consistently with 
the discussion here.) One can write p^ = b\b j = &t®J 
and 



(6) 



so that to compute the one-particle density matrix pij 
(pij) one only needs to calculate 



i-i 3-1 
ft ; 



/3=1 7=1 



det rp*vp j 



where 



N L 

i^)=nE p ^4i°) 

K=l Q—l 



(7) 



(8) 



is the Slater determinant corresponding to the fermionic 
wave- function [L is the number of lattice sites), and 
(P a )L,JV+i, with a = is obtained using properties 
of Slater determinants and written as 



—P gK for q < a, K = 1, . . . , N 
Pqk, for q > a, k = 1, . . . ,N 
5 ag for k = N + 1 



(9) 



Once p^ is computed, the momentum distribution 
function can be determined using the Fourier transform 



Pij, 



(10) 



where a is the lattice constant. 



where K = 2n/a is the reciprocal lattice vector and n is a 
nonzero integer. The second and third terms in Eq. 
can be computed using the approach mentioned in the 
previous subsection, so we focus here on how to compute 
the first term 



(n k n k > 



l_ e lka{l - ])+lk ' a{l - m \b\bib\b m ), (12) 



L 2 



ijlra 



for which we extend the recipe for calculating the two- 
point correlations [25|, [26[ to obtain four-point correla- 
tions and hence the noise correlations. 

From the mapping between hard-core bosons and spins 
one gets the following expression for the four-point cor- 
relation function in terms of spin operators 

(b\bjb\b m ) = 5 3l {b\i m ) + {b]b]bjb m ) 

= 25 jl (ai*-) + (-!)*» {&+VT&+&-) (IS) 

where in the last step we have used Eq. ([6]). 

Next we note that last term in Eq. (fTB"]) can be rewrit- 
ten as 

= Sij5i m + (—l) Sii Si m Gji + (—l) Slm 5ijG m i 
+(-l) 5ii+Slm 5 im Gji + (-l) f « + *"» +tf «»Gf imil , (14) 



where G 



ijkl 



Note that all Gu can be 



obtained as described in the previous subsection. 

Using the Jordan- Wigner transformation in Eq. (01 , 
the four-point Green's function for the spin-1/2 system 
can be written as 



i-l 



3-1 



G 



ijkl 



(^in^ ; "Aii e "^ 



a=l 
fe-1 



9=1 
i-l 



/ II' 

7=1 



™ftL ft 



flUe-^^p), (15) 



5=1 



which using properties of Slater determinants, as de- 
scribed in Refs. [25|,[26j], can be computed as 



G 



ijkl 



det 



(16) 



where (P Q ^) l,n+2, with a{(3) = k, I, is given by 



B. Noise correlations 

In this work we are interested in the second-order 
correlations of hard-core boson systems in the quasi- 
momcntum space 28[. These noise correlations are de- 
fined as 

A kk > = {n k n k i) - {h k ){h k >) - (h k ) 8 k - k >, n K, (H) 



pap 



-PP K for q< a, k = 1, . . . , N + 1 
P$ K iov g>a,K = l,...,N + l (17) 
5 ae for k = N + 2 



and (P ,3 )l,jv+i is given by Eq. ([9]). This means that to 
determine each element of the four-point Green's func- 
tion we need to multiply a matrix of dimension (iV+2) x L 
by a matrix of dimension L x (N + 2) [an operation that 
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scales as (JV + 2) 2 L] and then compute the determinant 
of the resulting (JV + 2) x (JV + 2) matrix [an operation 
that scales as (JV + 2) 3 ]. Finally, to compute the full four- 
point Green's function, we need to calculate of the order 
of L A nonzero elements; that is, the total computation 
time scales as L 4 [A(N + 2) 2 L + B(N + 2) 3 ], with A and 
B being prefactors for matrix multiplications and matrix 
determinants, respectively. 

III. HOMOGENEOUS CASE 

In this section we study the scaling of noise correla- 
tions in homogeneous chains. We should stress that, for 
all hard-core-boson systems considered in the following, 
periodic boundary conditions are always implemented; 
that is, for the equivalent fermionic Hamiltonians, peri- 
odic or antiperiodic conditions are selected depending on 
the number of particles in the lattice. 

In Fig. [TJ we show a typical noise correlation pattern 
for a strongly interacting supcrfluid system. It was cal- 
culated in a periodic lattice with 200 sites at half-filling. 
There are three features in that pattern that are appar- 
ent. First, a very large peak appears at k = k' = 0, re- 
flecting the presence of quasicondensation in the system. 
Replicas of this peak also appear at integer multiples of 
the reciprocal lattice vector K. Second, a line of max- 
ima can be found for k = k' due to the usual bunching 
in bosonic systems. Finally, dips are seen along the lines 
k, and 0, k' , which are related to the quantum depletion 
in the system. These features have been discussed in de- 
tail by Mathey et at. [3^| in the more general context of 
Luttinger liquids, for which hard-core bosons correspond 
to a limiting case. 
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FIG. 1: Noise correlations as a function of k and k' for a 
homogeneous system with 100 hard-core bosons in 200 lattice 
sites. 

As a function of the density p = N/L, the evolution 



of the noise correlations along the line k, is depicted in 
Fig. [2] The dips around the k = 0, ±K peaks are more 
clearly seen in Fig. [2j As noted in Ref. [33|, we find that 
the maximum value of Aoo occurs for p > 0.5, making 
evident the breakdown of the particle-hole symmetry for 
this observable. 
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FIG. 2: Noise correlations for k' = as k and p are changed 
for a system with L = 200. 

In what follows, we will focus on the scaling of the Aoo 
for different densities. For the k = peak of the mo- 
mentum distribution function, it is well known that the 
power-law decay of the one-particle correlations results 
in a rifc = o ~ s/L scaling (la] . In Fig. [3l we show the 
scaling of Aoo for three different densities in our periodic 
systems. 




FIG. 3: Scaling of the noise correlations Aoo for three different 
densities, p =0.25, 0.5 and 0.75. The solid lines are numerical 
fits in the form of Eq. (|18[) and the results (see text) show a 
leading-order linear behavior. 

To numerically find the scaling with system size, we 
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assume that 

A Qa = aL x + bL y , (18) 

where x and y (y < x) describe the leading and sublcad- 
ing terms, respectively, and a and b are coefficients that, 
together with x and y, are determined by means of a 
numerical fit. The results obtained for those four fitting 
parameters are given in Table HI 

TABLE I: Fitting parameters for the homogeneous case. 





p = 0.25 


p = 0.5 


p = 0.75 


a 


0.17779(2) 


0.21(9) 


0.16(1) 


X 


1.00041(2) 


1.010(3) 


1.010(8) 


b 


-1.049(4) 


0.15(7) 


0.534(6) 


y 


-0.099(1) 


0.8(1) 


0.59(2) 



Table Q] shows that the leading term is essentially lin- 
ear in all cases, while the exponent of the power law of 
the subleading term does depend on the density and was 
found to be quite close to one around half-filling. This 
means that in ultra-cold gas experiments one would need 
to reach large systems sizes to be able to clearly observe 
the linear scaling of the noise correlation peaks around 
half-filling, while this scaling would be more easily ob- 
served far away from half-filling. 

IV. PERIOD-TWO SUPERLATTICE 

We now consider the case in which an additional lat- 
tice, with twice the periodicity of the original lattice, is 
added to the system (a superlattice) . In this case, the 
on-site potential in Hamiltonian ([1} has the form 

Vi = Vcos(m), (19) 

with V representing its strength. As discussed in Ref. [36|, 
the effect of a period- two superlattice is to open a 
gap of magnitude V in the energy spectrum, splitting the 
original band into two bands. As a result, besides the 
usual insulating phases at p = 0,1, the half- filled system 
in the ground state also exhibits insulating behavior. As 
opposed to the p = 0, 1 insulators, the p = 0.5 (Mott) 
insulator does exhibit nonzero density fluctuations and a 
finite correlation length. 

Figure [4] shows the noise correlation pattern for the 
p = 0.5 insulator with V = It. Broad peaks can be 
clearly seen along the line k = k' , and those are character- 
istic of the noise correlations in the fractional Mott phase. 
They contrast with the sharp peaks seen in the noise cor- 
relations of the superfluid regime studied in the previous 
section. The suppressed height of the peaks in Fig. @] 
is a signature of the destruction of the quasi-long-range 
coherence in the half-filled Mott system. At this critical 
filling, the power-law decay of the one-particle correla- 
tions present in the absence of the superlattice is substi- 
tuted by an exponential decay pij ~ cxp(— \i — for 




FIG. 4: Noise correlations for the fractional Mott phase in the 
half-filled system in the presence of period-two superlattice for 
L = 200. 

which the correlation length £ was found to be £ ~ Ijy 
for small values of V (V < t) and £ ~ 1/ y/V for large 
values of V (V > t) [33]. As long as the lattice sizes 
are sufficiently large (L>Q, the absence of quasi-long- 
range coherence should be clearly observed in the scaling 
of the noise correlations in those systems. 

In the presence of the superlattice potential, additional 
features emerge in the noise correlations for k = k' ±n/a. 
Those can actually be used to distinguish the fractional 
insulator from the integer Mott insulator as both have 
suppressed Aoo peaks but only the former has a structure 
in the noise correlations for k = k' ± n/a. 




FIG. 5: Noise correlations A^o as a function of k and the 
density for superlattice systems with (a) V = It, (b) V = 2t, 
(c) V = 3t, and (d) V = 4t. L = 200 in all cases. 

In Fig. [5J we present a unified view of the behavior of 
the noise correlations for different systems with a super- 
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lattice potential. There we plot Aj-o as a function of p 
and k for four different values of V. Figure [5] shows that 
as V increases from It to At, the intensity of the central 
peak decreases for all fillings. However, the- suppression is 
more dramatic around half-filling. The additional unique 
signature of the presence of a superlattice potential is the 
structure that can be found at ka = ±7r. It is usually a 
positive peak for densities below 0.5 and becomes a dip 
right after the density increases beyond the fractional 
filling insulating phase. This peak-to-dip transition was 
discussed in detail by Rey et at [U, [3g , where in the 
limit V — > oo one can show analytically that Aoo and 
A±no have different signs for N = L/2 + 1. 





FIG. 7: Scaling of Aoo for p = 0.25 and p = 0.75 in sys- 
tems with V = It and lattices with 200 sites. The top inset 
shows results for the same systems in the fractional (p = 0.5) 
Mott-insulating phase. The bottom inset shows scaling of the 
sublattice peak Azlo for p = 0.25 in the same systems. Solid 
lines are numerical fits to the results, exhibiting a leading 
linear scaling with L in all superfluid cases. 



Assuming the same scaling ansatz in Eq. (jl8p. but in 
the presence of the superlattice only used for the super- 
fluid phases, we obtain the values depicted in Table [TT1 for 
the fitting parameters 



FIG. 6: The sublattice peak A» ( o as a function of p for four 
different values of V in systems with 200 sites. The inset 
shows the same quantity for systems with 100 sites. 

The behavior of A no as a function of the density and 
for different values of V can be better seen in Fig. |6] 
Interestingly we find that, in addition to the peak to dip 
transition around p = 0.5, there are other dip to peak and 
peak to dip transitions for higher densities. Those are 
only apparent for sufficiently large system sizes (beyond 
the ones studied in Refs. [3J, |35(). The inset in Fig. [S] 
shows that for a smaller system size with only 100 sites 
Ajlo is always negative for p > 0.5. 

Now that the generic features of the noise correlations 
in a superlattice potential have been reviewed, we focus 
on the scaling of the peaks with system size. In the frac- 
tional insulating regime, one expects that the exponential 
decay of correlations should lead to a saturation of the 
noise correlation peaks. This is, indeed, what we find, 
and an example is depicted in the top inset in Fig. [7] for 
half-filled systems with V = t. 

In the superfluid phases, on the other hand, it has 
been shown that one-particle correlations decay with ex- 
actly the same power law as the homogeneous system 
[3?| . Hence, we expect to find the same leading order 
scaling of A o that was discussed for homogeneous sys- 
tems in the previous section. This result can be seen in 
the main panel of Fig. and it can also be seen for the 
A»o peak, for p = 0.25, in the bottom inset. 



TABLE II: Fitting parameters for the superlattice case. 





k = 





k — n/a 




p = 0.25 


p = 0.75 


p = 0.25 


a 


0.15928(6) 


0.090(1) 


0.00716(2) 


X 


0.99023(6) 


1.045(1) 


1.0133(3) 


b 


-5.2(3) 


0.5793(5) 


-0.03891(9) 


y 


-0.71(1) 


0.636(2) 


0.417(1) 



As expected, we find that the leading terms of the noise 
correlation peaks are also of order L for both Aoo and 
Ailo- Similar to the homogeneous case, the leading lin- 
ear scaling of those peaks is better seen at low densities 
where finite-size effects have been found to be smaller 
because the subleading term has a much slower scaling 
with system size than the leading term. 



V. DISORDERED SYSTEM 

The disordered case is simulated by a random on-site 
potential of the form 

Vi = Se u (20) 

where S represents the strength of disorder and {et} are 
a set of random numbers between -1 and 1 selected with 
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a uniform probability distribution. For our disorder cal- 
culations we usually average over between 128 and 256 
disorder realizations. 

For one-dimensional noninteracting fermionic systems, 
the presence of disorder is known to lead to Anderson 
localization. This is a phase in which correlations de- 
cay exponentially while the system remains compress- 
ible; that is, no gap is present in the energy spectrum. 
Since hard-core bosons can be mapped to noninteracting 
fermions, the same is known to be true for the former. 
We should note that despite the fact that the one-particle 
correlations of hard-core bosons are in general different 
from those of noninteracting fermions, they also decay 
exponentially. This is shown in Fig. [8j where we present 
p x (with x = \i-j\) for systems with different disorder 
strengths. One should note that the exponential decay 
always sets in beyond a certain distance, which decreases 
as the strength of the disorder increases; that is, small 
systems with weak disorder may behave as superfluids. 
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FIG. 8: The decay of one-particle density matrices in half- 
filled systems with L = 500, characterized by different disor- 
der strength. The disorder averaging is performed over 128 
realizations for all S cases. Solid lines depict exponential 
decay, except for the homogeneous (5 = 0) system, where a 
solid line depicts the known power law y/x. Note the log-linear 
scale. 

From the previous discussion and the results in Fig. [8] 
one expects that, for any given system size, the height 
of the noise correlation peaks should decrease with in- 
creasing disorder strength and the peaks should become 
broader. This can be seen in Fig. [9j where we show 
the noise correlations A^i for four different disordered 
strengths in systems with 100 sites. The pattern for the 
S = 1 case resembles that of a homogeneous supcrfluid 
system (Fig. [T]), while for larger values of S they dis- 
play more similarities with the fractional Mott-insulating 
phase in the half-filled superlattice systems, with a clear 
broadening of the peaks at k = k'. (Of course, no addi- 
tional feature appear for k = k! ± ir/a in the disordered 
case.) 

A comparison between the cross-sectional view (for 
k' = 0) of the noise correlations in all three phases 




FIG. 9: Disorder-averaged noise correlations as a function of 
k and k' for systems with different disorder strength S =(a) 
It, (b) 2t, (c) 3t, and (d) At. N — 50 and L = 100 for all cases 
and the average was performed over 128 disorder realizations. 

discussed previously, namely, the superfluid, fractional 
Mott, and glassy phases, is shown in Fig. [TO] This com- 
parison makes evident (i) the suppression of the Aoo 
peak in the fractional Mott and glassy phases, (ii) the 
fact that the two insulating phases can be distinguished 
by the superlattice-induced features at k = k' ± it, and 
(iii) that the disordered Anderson-glass and the super- 
fluid phase exhibit the same satellite dips accompanying 
the Aoo peaks, while the dips vanish rather quickly in 
the fractional Mott phase. 
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FIG. 10: Noise correlations with fixed k' = for three half- 
filled systems with L = 200. The superfluid phase, the frac- 
tional Mott phase, and the Anderson-glass phase are associ- 
ated with the homogeneous, period-two superlattice (V — 2t), 
and disordered (5 = 2t) cases, respectively. The average is 
performed over 256 disorder realizations. 

Similar to the behavior of the fractional Mott phase, 
one also expects that as the system size increases for 
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any nonzero value of the disorder strength, the Aoo will 
saturate to a size-independent value that will only be a 
function of the density and the disorder strength. This 
behavior is shown in Fig. QT] for three different values 
of the disorder strength and for systems with up to 200 
sites. 
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FIG. 11: Scaling of the noise correlations Aoo for three dif- 
ferent values of 5 in half-filled systems. 

Finally, we study how the Aoo peak in the noise corre- 
lations behaves as a function of the disorder strength for 
a fixed size of the lattice. Since for S — we have already 
shown that such a peak diverges with system size, in the 
following we analyze how Aoo decreases as the disorder 
strength increases. 

In Fig.[12l we show Aoo as a function of 5 for two differ- 
ent system sizes. Three different regimes can be clearly 
identified, (i) For small values of 5, Aoo approximately 
stays constant with the increase of S, which can be un- 
derstood to be a consequence of a correlation length that 
exceeds the system size. As seen in FiglTSl that region 
decreases as the system size increases, (ii) As S increases 
even further, a power- law decay develops in Aoo, and the 
region over which such a power law can be seen increases 
with system size as regime (i) is suppressed. In our fits, 
we find the power law Aoo ~ to have an exponent 
7 ~ 1.78(2), but it is still influenced by some finite-size 
effects. In order to gain further understanding of the 
power-law decay of the height of this noise peak in insu- 
lating phases, we have studied the behavior of A o vs V in 
the fractional Mott phase in a superlattice, for which we 
can study larger systems sizes. We find that Aoo ~ V ~ 7 
with an exponent of 0.874(5), which is different from the 
one for the disordered system. These results clearly show 
that the power-law decay of Aoo as one enters an insu- 
lating phase depends on the perturbation creating the 
insulator, i.e., it is not universal, (hi) Finally, for very 
strong disorder, Aoo saturates to a nonzero value. This 
asymptotic behavior is found to agree with the analytical 
value in the 5 — > oo limit, computed using 

A 00 = p(p+f), (21) 



which was derived by Rey et al. [34|. Equation (f2"Tj) 
shows that A o only depends on the density and also 
makes explicit the absence of particle-hole symmetry for 
this observable in hard-core-boson systems. This third 
regime is robust against the disorder variance, something 
that follows from the fact that the correlation length is 
of the order of or smaller than the lattice spacing a. 




FIG. 12: Noise correlation peaks Aoo as a function of disorder 
strength in two half-filled disordered systems with L — 100 
and 200, respectively. The solid line shows a power-law fit 
Aoo ~ S with 7 = 1.78(2) in the range from 8 — t to 
3t for L — 200. The top inset shows Aoo as a function of 
V in three half-filled period-two superlattices. The solid line 
depicts a power law with an exponent 7 = 0.874(5). The 
bottom inset shows the asymptotic behavior for large values 
of 5; the dotted line marks the analytical result in the limit 
of infinite disorder. 



VI. CONCLUSIONS 

We have implemented an exact approach to numeri- 
cally compute the noise correlations for hard-core boson 
in one-dimensional lattices. For that purpose, we have 
extended to four-point correlations the recipe for calcu- 
lating two-point correlations introduced in Refs. [25l. l26j . 
Our approach has a polynomial time scaling that is more 
efficient than the straightforward application of Wicks 
theorem, and can be easily extended to study noise cor- 
relations in nonequilibrium systems. 

We have applied this approach to study the scaling of 
noise correlations in three different phases that appear in 
homogeneous systems and in the presence of two differ- 
ent background potentials. We have shown that in the 
supcrfluid phase, the noise correlation peaks Aoo exhibit 
a leading linear behavior ~ L, independent of the density 
and of the presence of a superlattice potential. The sub- 
leading term was found to be strongly dependent on the 
latter two. On the other hand, the fractional Mott and 
Anderson-glass insulating phases exhibit an asymptotic 
value, which is independent of system size and only de- 
pends on the density and the strength of the potentials 
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creating such phases. This behavior was expected and 
manifests the absence of quasi-long-range order in these 
two phases. 

In the period-two superlattice, we have also found vari- 
ous peak-to-dip and dip-to-peak transitions that were not 
observed in previous studies with smaller system sizes, 
something that demonstrates the importance of finite- 
size effects in the noise correlations and the need for ap- 
proaches that allow one to study very large system sizes. 

Finally, we have shown that in the disordered system 
(fractional Mott phase), the decrease of the A o peak 
with increasing disorder strength (superlattice strength) 
exhibits a region with a power-law decay Aqo ~ S(V) 1 , 



with a nonunivcrsal value of exponent 7 that depends on 
the kind of perturbation creating the insulator. 
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